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Abstract 

Rhagoletis cerasi (Diptera: Tephritidae) is a major pest of sweet and sour cher- 
ries in Europe and parts of Asia. Despite its economic significance, there is a 
lack of studies on the genetic structure of R. cerasi populations. Elucidating the 
genetic structure of insects of economic importance is crucial for developing 
phenological-predictive models and environmental friendly control methods. 
All natural populations of R. cerasi have been found to harbor the endosymbi- 
ont Wolbachia pipientis, which widely affects multiple biological traits contrib- 
uting to the evolution of its hosts, and has been suggested as a tool for the 
biological control of insect pests and disease vectors. In the current study, the 
analysis of 18 R. cerasi populations collected in Greece, Germany, and Russia 
using 13 microsatellite markers revealed structuring of R. cerasi natural popula- 
tions, even at close geographic range. We also analyzed the Wolbachia infection 
status of these populations using 16S rRNA-, MLST- and wsp-based approaches. 
All 244 individuals screened were positive for Wolbachia. Our results suggest 
the fixation of the wCerl strain in Greece while wCer2, wCer4, wCer5, and 
probably other uncharacterized strains were also detected in multiply infected 
individuals. The role of Wolbachia and its potential extended phenotypes needs 
a thorough investigation in R. cerasi. Our data suggest an involvement of this 
symbiont in the observed restriction in the gene flow in addition to a number 
of different ecological factors. 
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Introduction 

The European cherry fruit fly, Rhagoletis cerasi L. (Dip- 
tera: Tephritidae) (Fig. 1), is a pest of major agricultural 
importance infesting mainly sweet (Prunus avium) and 
sour (P. cerasus) cherries as well as honeysuckle (mainly 
Lonicera tatarica and L. xylosteum) and wild growing 
prunus. Rhagoletis cerasi has dispersed and become estab- 
lished in almost all European countries, spreading from 
the original habitat, the Caucasian area of western Asia 
(Fimiani 1989; White and Elson-Harris 1992). Rhagoletis 
cerasi follows a univoltine life cycle with obligatory pupa 
dormancy that lasts over 9 months. Adults emerge in 
spring and oviposit usually a single egg on ripe or ripen- 
ing fruits (Fletcher 1989; Daniel and Grunder 2012). 
Larvae feed on mesocarp, destroying fruits and causing 
considerable economic loss exacerbated by secondary fun- 
gal and bacterial infections (Fimiani 1989). Because of the 
economic importance, the control of the European cherry 
fruit fly received much attention in the framework of the 
sterile insect technique (SIT) in the early 1970s (Boiler 
et al. 1971, 1977; Boiler and Katsoyannos 1977). 

A number of studies have focused on the ecology of 
R. cerasi (Boiler and Prokopy 1976; Boiler et al. 1998; 
Kovanci and Kovanci 2006; Papanastasiou et al. 2011; 
Moraiti et al. 2012a,b) or the development of environ- 
mental friendly control methods, such as mass trapping, 
application of the oviposition deterring pheromone, bait 
sprays using plant-derived insecticides and biological 
control (Katsoyannos et al. 2000; Daniel and Grunder 2012). 
However, very little is known regarding the biology of 
this pest at the molecular, genetic, and population levels. 
The cytogenetics (including the development of polytene 
chromosome maps) have been recently studied (Kounati- 
dis et al. 2008), though the genetic structure of the Euro- 
pean cherry fruit fly populations has received much less 
attention. Allozyme polymorphisms (Schwarz et al. 2003) 
raise the question of host race formation in this species, 
due to host shift, and a small-scale analysis based on 
microsatellite markers (Augustinos et al. 2011a) points 
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toward elevated genetic differences even in a short 
geographical range. 

In fruit flies (other than Rhagoletis species), microsatel- 
lites have been used to shed light on: recent invasion 
events and dispersion routes of several invasive species 
(Meixner et al. 2002; Bonizzoni et al. 2004; Augustinos 
et al. 2005; Nardi et al. 2005; Aketarawong et al. 2007; 
Khamis et al. 2009; Zygouridis et al. 2009; Virgilio et al. 
2010); the degree of genetic structuring in natural popula- 
tions that can be associated with incipient speciation 
(Michel et al. 2007); and pattern of remating in nature 
(Bonizzoni et al. 2002; Song et al. 2007). Lately, microsat- 
ellite markers have been also used to address phylogenetic 
issues in the Bactrocera musae complex of species (Drew 
et al. 2011). 

Microsatellite markers have been recently developed for 
R. pomonella (Velez et al. 2006), R. indifferens (Maxwell 
et al. 2009), and R. cerasi (Arthofer et al. 2009a) and 
used to: (a) address issues regarding speciation events in 
R. pomonella (Michel et al. 2007, 2010; Forbes et al. 
2009), and (b) analyze native and introduced populations 
of R. completa (Chen et al. 2010). The recent develop- 
ment of 13 microsatellite markers for R. cerasi (Arthofer 
et al. 2009a), along with the evaluation of cross-amplified 
microsatellite markers from other Rhagoletis species 
(Augustinos et al. 2011a), provide a solid basis for explor- 
ing the population genetics of this species. 

Wolbachia pipientis (Wolbachia for the purpose of this 
manuscript) is perhaps the most widespread and abun- 
dant symbiont among insect species. Extensive genetic 
diversity has been recorded in the Wolbachia group with 
the different strains assigned in 13 supergroups, namely A 
to F and H to N (Augustinos et al. 2011b). The genotyp- 
ing of the Wolbachia strains is based on either single gene 
(such as wsp, ftsZ, gltA, groEL, dnaA) or Multi Locus 
Sequence Typing (MLST) approaches (Casiraghi et al. 
2005; Baldo et al. 2006; Paraskevopoulos et al. 2006). 

The widespread distribution of Wolbachia, its ability to 
manipulate reproductive properties, as well as other 
important physiological functions of its arthropod hosts, 



Figure 1. Rhagoletis cerasi fly: (a) third instar 
larvae of R. cerasi infesting a cherry fruit and 
(b) R. cerasi adult (male) on food droplets. 
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have attracted the attention of investigators regarding the 
role of this symbiont in host biology, ecology, and evolu- 
tion (Serbus et al. 2008; Werren et al. 2008; Saridaki and 
Bourtzis 2010). In addition, Wolbachia symbiosis is 
currently being explored and harnessed toward population 
control of insect pests and disease vectors (Zabalou et al. 
2004, 2009; Brelsfoard and Dobson 2009; Apostolaki et al. 
2011; Atyame et al. 2011; Walker et al. 2011). 

Revealing the close association between R. cerasi and 
Wolbachia, Riegler and Stauffer shed important light on 
the reproductive incompatibility reported in the mid- 
1970s between northwest and southeast European popula- 
tions of this pest (Boiler et al. 1976; Matolin 1976; Riegler 
and Stauffer 2002). So far, five Wolbachia strains of 
R. cerasi have been described (wCerl to wCer5). The 
identification of the first two strains, wCerl and wCer2, 
was based mainly on wsp gene divergence. Wild _R. cerasi 
populations were found to be infected either with wCerl, 
or super infected with wCerl and wCer2. Moreover, 
wCer2 was proposed to be the causal factor for the repro- 
ductive incompatibility of the aforementioned popula- 
tions, now known as Wolbachia-induced cytoplasmic 
incompatibility (CI) (Riegler and Stauffer 2002). Based on 
a polymerase chain reaction (PCR) screen with Wolbachi- 
a-specific primers for 16S rRNA and wsp, Kounatidis 
et al. (2008) reported that the four Greek populations 
tested were also 100% infected. The Wolbachia strains 
wCer3, wCer4, and wCer5 were only recently isolated 
from R. cerasi populations (Arthofer et al. 2009b, 2011). 
Interestingly, wCer4 was identified only after Ceratitis 
capitata transinfection with Wolbachia, using _R. cerasi as 
donor (Zabalou et al. 2004). Arthofer et al. (2011) pro- 
posed Allelic Intersection Analysis as a novel tool for the 
MLST assignment of multiply infected hosts and 
described the MLST profile of R. cerasi Wolbachia strains 
wCerl to wCer5 (except wCer3), indicating that wCerl, 
wCer2, and wCer4 belong to Wolbachia supergroup A, 
while wCer5 belongs to supergroup B. The Wolbachia 
strain wCer3 is supposed to be a recombinant A/B strain, 
according to the available wsp gene sequence data (Artho- 
fer et al. 2009b). 

Recent studies showed that R. pomonella (Schuler et al. 
2011) and R. cingulata (Drosopoulou et al. 2011) are 
infected by Wolbachia, while _R. completa seems to be free 
of Wolbachia (Drosopoulou et al. 2010). Rhagoletis pomo- 
nella is infected with the wPom strain, which seems to be 
very similar (or identical) to wCer2, while other sequence 
variants for different MLST genes (not attributed to 
specific strains) have also been reported (Schuler et al. 
2011). Rhagoletis cingulata is infected with at least two 
strains very similar (or identical) to wCerl and wCer2 
(Schuler et al. 2009; Drosopoulou et al. 2011), named 
wCinl and 2 (Schuler et al. 2009). 
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Wolbachia strains derived from R. cerasi have been 
recently used for the development of an alternative and 
environment-friendly strategy to control two major agri- 
cultural pests: the Mediterranean fruit fly, C. capitata 
(Zabalou et al. 2004, 2009) and the olive fly, Bactrocera 
oleae (Apostolaki et al. 2011). This technique, the Incom- 
patible Insect Technique (IIT), is in principle analogous 
to the SIT. The main difference between them is the ster- 
ilization factor; the IIT is based on Wolbachia-induced CI 
while the SIT is based on irradiation (Bourtzis and Rob- 
inson 2006). IIT holds great potential for the control of 
agricultural pests and disease vectors (Zabalou et al. 2004; 
Brelsfoard and Dobson 2009; Saridaki and Bourtzis 2010; 
Apostolaki et al. 2011; Atyame et al. 2011). 

Because of its wide geographical distribution, univoltine 
life history, patchy distribution of its hosts, low adult dis- 
persion abilities, and infection with different Wolbachia 
strains (Phillips and Dirks 1933; Jones and Wallace 1955; 
Boiler and Prokopy 1976; Fletcher 1989; Kneifl et al. 
1997), R. cerasi bears several properties of becoming an 
important model for addressing issues regarding effect of 
ecological factors and symbionts on the differentiation of 
genetic, demographic, and behavioral traits (Augustinos 
et al. 2011a; Moraiti et al. 2012b). 

Materials and Methods 

Collection of fly samples and DNA 
extraction 

A total of 465 adults from 18 different R. cerasi popual- 
tions were genotyped using 13 microsatellite markers. 
Fifteen of the populations came from Greece (12 main- 
land and 3 island) and additionally three from Germany 
(two populations) and Russia. In all cases, field-infested 
sweet cherries were sampled and kept in the laboratory 
until the collection of pupae. Adults emerging from 
pupae were stored at — 80°C in 95% ethanol. Genomic 
DNA was extracted as described in Augustinos et al. 
(2011a). Collection sites and the number of flies used in 
the study are shown in Figure 2 and Table 1. 

Microsatellite genotyping and data analysis 

PCRs, genotyping, and specific allele sequencing of micro- 
satellite markers were performed as described in Augusti- 
nos et al. (2011a). PCR conditions for markers Rce76-1 
and 83-44 were similar to those described in Arthofer 
et al. (2009a). 

Genetic variability was measured as the mean number 
of alleles per locus (« a ), effective number of alleles (« e ), 
observed (H 0 ) and expected heterozygosity (H e ), using 
POPGENE, version 1.31 (Yeh et al. 1999). Allelic diversity 
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Figure 2. Rhagoletis cerasi sampling areas and GeneClass analysis. Sampling sites: 1, Agia Larissa; 2, Kallipefki Larissa; 3, Kato Lechonia; 4, 
Kamari; 5, Karditsa; 6, Pertouli; 7, Konitsa; 8, DafnM; 9, Dafni2; 10, Kastoria; 11, Salonica; 12, Kernitsa; 13, Chania; 14, Chios; 15, Lesvos; 16, 
Steckelberg; 17, Dossenheim; 18, Krasnodar. Samples were grouped for GeneClass analysis according to geographic origin. 



was determined after correction for sample size, using 
FSTAT (Goudet 2001). Deviations from Hardy-Weinberg 
Equilibrium (HWE) were tested with the G 2 likelihood 
ratio test in POPGENE. Genotypic disequilibrium was 
tested with Genepop (Raymond and Rousset 1995), using 
Fisher's exact test, for all pairs of loci in all populations 
and across them. Presence of null alleles was tested with 
ML-Null software. Hardy-Weinberg tests for heterozygote 
deficiency were performed (Rousset and Raymond 1995), 
after performing 10,000 randomization tests using Monte- 
Carlo randomization as described by Guo and Thomson 
(1992). Genetic distances were measured according to Nei 
(1972), using POPGENE (Yeh et al. 1999). Population 
differentiation was estimated using the FSTAT software, 
using a significance level of 0.05, adjusted for multiple 
comparisons. Software GENALEX 6.5 (Peakall and 
Smouse 2006) was used for the estimation of the pairwise 
population PhiPT values. This PhiPT matrix was used in 
Genalex 6.5 to perform Principal Components Analysis 
(PCA). PHYLIP 3.6c (Felsenstein 1994) was used for the 
construction of an unrooted UPGMA dendrogram after 
1000 bootstrap resamples, using allele frequencies, as 
described in Augustinos et al. (2005). STRUCTURE v2.2 
software was used to determine the number of possible 



genetic clusters (Pritchard et al. 2000). We used all four 
different models implemented in this program, assuming 
the presence of admixture and without admixture, for 
independent or correlated allele frequencies, with a burn- 
in period of 50,000 and 500,000 Marcov Chain Monte 
Carlo (MCMC) repetitions after the initial burn-in. We 
tested for K = 1 to K = 12 (where K stands for the 
assumed number of populations). We also used the mod- 
ification described by Evanno et al. (2005) to predict the 
true number of clusters in complex situations more accu- 
rately. For this purpose, we used the "no admixture" 
model with correlated frequencies, with a burn-in period 
of 50,000 and 500,000 MCMC repetitions after the initial 
burn-in. We ran 10 repetitions of the above model, 
assuming K = 1 to K = 12. We chose the specific model 
because, as it is suggested by the authors of the program, 
it can better resolve "complicated" or "weak" cases of 
structuring. Software GeneClass 2.0 (Piry et al. 2004) was 
used to perform population assignment and exclusion test 
and to calculate the probability of origin for each individ- 
ual and each population, under a Bayesian model using 
the Rannala and Mountain (1997) criterion with a 0.05 
threshold. Isolation by distance was tested using the 
"Mantelize it" option of the FSTAT software (Goudet 
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2001) for plotting genetic distances against geographic 
distances, and Pearson's r coefficient was measured using 
Excel. Finally, the BOTTLENECK software was used 
(Cornuet and Luikart 1996) to detect any recent bottle- 
neck phenomena. This software runs under the assump- 
tion that recent bottlenecks lead to a shift from the L- 
shaped distribution of allele frequencies along with a 
quicker loss of rare alleles than heterozygosity. However, 
the low number of microsatellite loci used, especially if 
not highly polymorphic, can lead in reduced "ability" to 
detect bottleneck phenomena. In the frame of this study, 
the Mode shift and the Wilcoxon sign rank tests imple- 
mented in this software were used. 

Wolbachia genotyping and data analysis 

Wolbachia infections from six individual flies belonging to 
separate genetic clusters (according to microsatellite analy- 
sis) and/or of different geographic origin were genotyped by 
sequencing of the 16S rRNA, wsp, and MLST genes. A 16S 
rRNA gene fragment (438 base pairs) was amplified with the 
Wolbachia-specific primers wspecF and wspecR (Werren and 
Windsor 2000). Gene fragments of wsp and the MLST genes 
(gatB, coxA, hcpA, fbpA, and ftsZ) were amplified using the 
respective primers and conditions reported previously 
(Casiraghi et al. 2005; Baldo et al. 2006; Ros et al. 2009). 

Cloning and sequencing 

After detecting multiple peaks in preliminary sequencing 
reactions, PCR products of all genes were cloned, sequenced, 
and assembled as described in Augustinos et al. (2011b). 
Three to six colonies were directly subjected to PCR, using 
the primers T7 and SP6 and double-strand sequenced. 

Phylogenetic analysis 

All Wolbachia gene sequences generated were aligned 
using the ClustalW algorithm (Thompson et al. 1994), 
implemented in MEGA (Tamura et al. 2011). Sequences 
obtained from GenBank representing all currently known 
supergroups of Wolbachia were included in the analysis. 
Phylogenetic reconstruction was performed using the 
Maximum Likelihood (ML) method, after 500 bootstrap 
resamples of the original data. 

Recombination analysis 

Sequences for all genes (MLST, wsp, 16S rRNA) were 
tested for recombination events using the RDP3 package 
(Martin et al. 2010). All the available modules (RDP3, 
BOOTSCAN, GENECONV, MAXCHI, CHIMAERA, 
SISCAN, 3SEQ, VisRD, LARD, and LDHAT) were 
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implemented under the default options. As the authors of 
the software state, the "default options" are sufficient to 
detect recombination events for most of the datasets even 
without the need of reference nonrecombinant sequences. 

Screening for wCer specific strains 

Ten to twenty flies per population were screened for the 
presence of wCerl to wCer5 Wolbachia strains using a 
wsp gene-based PCR assay as previously reported (Artho- 
fer et al. 2009b). 

Nucleotide sequence accession numbers 

All microsatellite and Wolbachia gene sequences generated 
in this study have been deposited into GenBank under 
accession numbers JX870457-JX870517 and as well as to 
the Wolbachia MLST database. 

Results 

Microsatellite marker variability and HWE 

The markers used have on average 6.23 n a with the mean 
w e being 2.16 (Table 2). The mean observed heterozygos- 
ity was lower than the expected (0.39 vs. 0.47). Out of 
the 234 locus-sample tests performed for HWE, 25 
showed deviation according to the G 2 criterion at a sig- 
nificance level of 0.05. Most of these deviations are con- 
centrated in marker Rpll (six) and can be partially 
attributed to the presence of null alleles, as this locus 
exhibits an excess of different classes of homozygotes. 
Analysis performed with the ML-Null software points to 
the presence of null alleles for this locus in some of the 
populations. No evidence for extensive presence of null 
alleles for the other loci was found. 

Population polymorphism and HWE 

Populations studied exhibit a varying degree of polymor- 
phism with 1.9 to 2.51 n e per locus (Dafnil-GR and 
Krasnodar-RU, respectively). According to all measures, 
Agia Larissa-GR and Krasnodar-RU are the most polymor- 
phic populations while Kernitsa-GR and the two neighbor- 
ing islands of Chios-GR and Lesvos-GR present the lowest 
variability. Sporadic HWE deviations (up to three loci) 
were observed in few of the populations studied (Table 1). 

Genetic distances and population 
differentiation 

Pairwise genetic distances (Nei 1972) among populations 
were determined ranging from 0.0140 to 0.1630 (Table 3). 
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Table 2. Genetic variability of microsatellite markers. 
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Marker 


Reference 


Size variation 


W 


"a 


n e 


AR 


Ho 




HWE 


RcMic(Boms3b) 


Augustinos et al. (2008) 


94-136 


462 


4 


1.34 


3.33 


0.24 


0.26 


2/18 


Rc(Rp1) 


Velez et al. (2006) 
Augustinos et al. (201 1) 


486^196 


433 


3 


1.88 


2.37 


0.26 


0.47 


2/18 


RcMic(Rp2) 


Velez et al. (2006) 


202-220 


456 


9 


2.96 


5.40 


0.63 


0.66 


1/18 


RcMic(Rpll) 


Velez et al. (2006) 


289-299 


416 


5 


2.60 


3.69 


0.40 


0.62 


6/18 


RcMic(Rp12) 


Velez et al. (2006) 


190-204 


463 


7 


2.47 


3.92 


0.50 


0.60 


4/18 


RcMic(Rp15) 


Velez et al. (2006) 


234-252 


451 


6 


1.90 


3.29 


0.43 


0.48 


0/18 


RcMic(Ri48) 


Maxwell et al. (2009) 


128-156 


454 


8 


1.50 


4.38 


0.22 


0.33 


1/18 


RcMic(Ri67) 


Maxwell et al. (2009) 


180-196 


445 


6 


1.42 


3.41 


0.27 


0.30 


0/18 


RcMic(Ri83) 


Maxwell et al. (2009) 


166-170 


461 


2 


1.46 


2.00 


0.30 


0.32 


2/18 


RcMic(Ri93) 


Maxwell et al. (2009) 


185-205 


449 


9 


3.70 


4.95 


0.57 


0.73 


3/18 


RcMic(Ri111) 
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3.54 


5.1 
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0/18 


Rce83-44 


Arthofer et al. (2009a) 


189-213 


459 


8 


2.00 


4.57 


0.40 


0.50 


0/18 


Mean 






452.5 


6.23 


2.16 


2.371 


0.47 


0.39 





W, sample size; n a , number of actual alleles; n e , number of effective alleles; AR, allelic richness; H 0 , observed heterozygosity; H e , expected hetero- 
zygosity; HWE, number of samples in which locus deviates from equilibrium, according to G 2 criterion, at a significance level of 0.05. 



fsT (fixation index) values were determined and tested at a 
significance level of 0.05 (Table 4). There are no significant 
differences among populations from Thessaly, except the 
one observed between Kamari Pilio-GR/Pertouli-GR. The 
four populations from Macedonia are also quite homoge- 
neous, with the exception of the difference between Kasto- 
ria and Dafni2. The single population from Epirus 
(Konitsa-GR) shows differentiation with some Thessaly/ 
Macedonia populations and at the same time, it is not sig- 
nificantly differentiated from the Dossenheim population. 
The populations from Peloponnesus (Kernitsa-GR), Crete 
(Chania-GR), Aegean islands (Lesvos-GR and Chios-GR), 
and Russia (Krasnodar-RU) are all differentiated from all 
previously reported and among them. The German popu- 
lations seem to form a homogeneous cluster. 

Analysis of the structure of R. cerasi natural 
populations 

To further unravel the genetic structure of R. cerasi natu- 
ral populations, different programs and assumptions were 
used. A UPGMA dendrogram was constructed in Phylip 
3.6c (Fig. 3). Bayesian clustering analysis was performed 
using STRUCTURE 2.2 (Fig. 4), along with the modifica- 
tion in Evanno et al. (2005), which is supposed to give a 
better resolution of the true number of clusters in complex 
cases of population structuring. PCA was also performed 
using Genalex 6.5 (Fig. 5). Although these approaches do 
not produce exactly the same results, the main conclusion 
is that there is some structuring with the populations from 
Krasnodar-RUS and Aegean Sea (Chios-GR and Lesvos- 
GR) to be strongly differentiated from all others. Differen- 
tiation is also evident for Kernitsa-GR (Peloponnesus) and 



the two German populations. Although there is not strong 
statistical support, two additional loose clusters may be 
present which include populations from Thessaly-GR and 
Macedonia/Epirus-GR. 

We used software GeneClass v2.0 to assign individuals to 
respective groups. Eight groups, based on geographic crite- 
ria, were formed: Thessaly-GR, Epirus-GR, Macedonia-GR, 
Peloponnesus-GR, Crete-GR, Aegean Islands-GR, Ger- 
many, and Russia. Some populations are well separated; the 
vast majority of their individuals are correctly assigned to 
the respective groups (Fig. 2). This is true for populations 
from Kernitsa-Peloponnesus (29/30 correctly assigned), the 
Aegean islands (56/60), Russia (28/32), and Stecklenberg 
(28/32). There was a high number of false assigned individ- 
uals for the remaining populations (ranging from 26% to 
more than 50% per population), especially for Macedonia- 
GR, Thessaly-GR, and Dossenheim-GER. 

Isolation by distance 

The hypothesis whether the observed genetic distances 
can be correlated with geographic distances between sam- 
pling locations was tested in pairwise comparison of f ST 
values and the logarithm of the geographic distances. 
There was only a weak correlation (r = 0,33, P < 0.1) 
indicating that other factors may also play a significant 
role in the development of the observed gene flow restric- 
tion between R. cerasi populations. 

Detection of bottlenecks 

The program BOTTLENECK was used for the detection 
of recent bottlenecks. Results have been handled with care 
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Figure 3. Unrooted UPGMA dendrogram created in Phyllip 3.6c, based on allele frequency differences, after 1000 bootstrap resamples. Only 
bootstrap values higher than BOO are shown. 




as only few of the populations consist of at least 30 indi- 
viduals and the loci used in this analysis are <20 and, in 
addition, they are not highly polymorphic. Given these 
limitations, no evidence for such phenomena was found 
with the exception of the Kernitsa population where a 
shifted mode of allele frequencies from the L- shaped 
distribution was observed. 



Wolbachia infections in R. cerasi 
populations 

Using a 16S rRNA gene-based PCR screen, all 244 individu- 
als originating from the 18 different R. cerasi populations 
were found infected with Wolbachia. To investigate whether 
there is a correlation between genetic structuring of 
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Figure 5. Principal Components Analysis 
(PCA), using the PhiPT distance matrix created 
in Genalex 6.5, based on the allele frequency 
differences among Rhagoletis cerasi samples. 
1, Agia Larissa; 2, Kallipefki Larissa; 3, Kato 
Lechonia; 4, Kamari; 5, Karditsa; 6, Pertouli; 7, 
Konitsa; 8, Dafnil; 9, Dafni2; 10, Kastoria; 11, 
Salonica; 12, Kernitsa; 13, Chania; 14, Chios; 
15, Lesvos; 16, Steckelberg; 17, Dossenheim; 
18, Krasnodar. 
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R. cerasi natural populations and the presence of different 
Wolbachia strains, we analyzed six individual flies from differ- 
ent groups as determined by microsatellite analysis and/or 
based on geographic criteria. These flies derived from the 
populations of Agia Larissa (Thessaly-GR), Chania (Crete- 
GR), Chios (Aegean Islands-GR), Kernitsa (Peloponnesus- 
GR), Krasnodar (Russia), and Stecklenberg (Germany). 

PCR products were obtained from these six flies for the 
16S rRNA, wsp as well as the five MLST genes (coxA, 
fbpA, ftsZ, gatS, and hcpA). All PCRs produced the 
expected amplicon and these products were cloned. Two 
to six clones were sequenced per product and were found 
to be Wolbachia specific: sixteen 16S rRNA, twenty-one 
coxA, twenty fbpA, twenty- two ftsZ, twenty gatS, 15 hep A, 
and eighteen wsp clones were sequenced revealing three to 
five alleles per locus (Table 5). The major findings of the 
sequencing analysis (Table 5) can be summarized as 
follows: (a) the majority of the 16S rRNA sequences (14/16) 
form a cluster in supergroup A (Figure S2). There are 
two exceptions: first, one sequence derived from the fly 
from Chania clusters in another clade, but still with 
sequences belonging to supergroup A and second, one 
sequence derived from the fly from Stecklenberg clusters 
with sequences belonging to supergroup B; (b) the 21 
coxA sequences revealed the presence of four alleles (Fig- 
ure S3). Most of the sequences (16) correspond to allele 
coxA-84, which has been reported in wCerl and wCer4 
strains, and were found in Greek populations (Agia, 
Chios, Kernitsa) as well as in Krasnodar and in Stecklen- 
berg. The allele coxA-1, which is present in wCer2, was 
detected in Stecklenberg too. Interestingly, two new alleles 
were detected in the individual from Chania-GR. These 
new alleles are named coxA-188 and coxA-189. CoxA-188 
differs only in 1 bp from allele coxA-130 and coxA-189 
differs in 1 bp from allele coxA-58. However, both cluster 



near the allele coxA-84, which is common in wCerl and 
VfCer4 strains. None of our coxA sequences group with 
wCer5, a Wolbachia strain belonging to supergroup B; (c) 
the 20 fbpA sequences suggested the presence of four 
alleles (fbpA-1, -4, -79, and -160) which correspond to 
alleles characterized so far to four different wCer strains. 
Most of the clones (15) correspond to allele fbpA-160, 
present in strain wCerl, and originate from Greek popu- 
lations (Agia, Chania, Chios, Kernitsa) as well as from 
Stecklenberg and Krasnodar. Two clones correspond to 
fbpA-1 of strain wCer2 and were detected in Stecklenberg, 
two clones correspond to fbpA-79, which is assigned to 
wCer4, and were present in Chania and in Kernitsa while 
one clone (from Chios-GR) corresponds to fbpA-4, which 
is found in wCer5, the only strain belonging to super- 
group B; (d) the 22 ftsZ sequences revealed the presence 
of four alleles (ftsZ-3, -22, -70, and -79) described for 
wCer strains. All sequences cluster with wCer strains 
belonging to supergroup A, with the exception of the sole 
sequence from Stecklenberg, which clusters with wCer5 of 
supergroup B (ftsZ-22). Thirteen of the clones correspond 
to the allele ftsZ-79 assigned to wCerl and originate from 
Greek populations (Agia, Chios, Kernitsa) and from 
Krasnodar. Six clones correspond to ftsZ-70 of wCer4 
(originating from Chania, Kernitsa, and Krasnodar) while 
two clones correspond to ftsZ-3 assigned to wCer2; (e) 
the 20 gatS sequences suggested the presence of three 
wCer alleles belonging to supergroup A (gatB-1, -8, and - 
53). Eleven of them correspond to allele gatS-8 of strain 
wCerl originating from Greek populations (Agia, Chios, 
Kernitsa) as well as from Stecklenberg and Krasnodar. 
Seven clones are identical to gatB-53 (wCer4) and were 
detected in Chania and Krasnodar while two clones are 
identical to the allele gatB-1, which is assigned to strain 
wCer2, and were found in Stecklenberg; (f) the 15 hcpA 
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Table 5. MLST and wsp alleles found in individual flies derived from Rhagoletis cerasi selected populations. 
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Gray shaded lines: MLST profile of wCer strains, as published by Arthofer et al. (201 1). 

In bold: new MLST alleles found in this study and deposited in the Wolbachia MLST database. 

MLST, multi locus sequence typing; N, number of colonies sequenced. 

In parentheses (): number of colonies harboring the respective allele. 

'in brackets []: supergroup to which each strain belongs. 



sequences revealed the presence of five wCer alleles 
(hcpA-40, -85, -103, -212, and -213). All of them cluster 
with wCer sequences of supergroup A, with the exception 
of two of the three clones from Chios, which are identical 
to allele hcpA-40, which is assigned to wCer5 (Figure S4). 
Eight clones correspond to allele hep A- 103 of strain 
wCerl originating from Greek populations (Agia, Chios, 
Kernitsa) and from Krasnodar. Only two sequences corre- 
spond to hcpA-85 of wCer4, both detected in Kernitsa. 
Interestingly, in the individual from Chania, two new 
alleles were present (-212, -213). HcpA-213 is very similar 
to hcpA-103 of strain wCerl (different in the first 3 bp), 
while hcpA-212 is very similar to hcpA-85 of strain wCer4 
(different in the first 3 bp also); and (g) the 18 wsp 
sequences suggested the presence of three known wCer 
alleles (wsp-113, -335, and -581). All of them cluster in 
supergroup A, with the exception of one of the five clones 
from Chios, which is identical to allele wsp-581, which is 
found in strain wCer5. Fifteen of the sequences were 
identical to wsp-335 of wCerl originating from Greek 
populations Agia, Kernitsa, Chania, Chios) and from 
Krasnodar, while two clones from Chania were identical 
to wsp-113 of strain wCer4. Taken together these data 
strongly suggest the presence of multiple Wolbachia infec- 
tions in all R. cerasi populations studied. 

The presence of multiple Wolbachia infections was also 
confirmed by a wsp gene-based diagnostic PCR assay pre- 
viously developed (Arthofer et al. 2009b). The primers 



used in this assay selectively amplify different alleles of 
the wsp gene that have been attributed to the different 
wCer strains. Ten to twenty individuals were screened per 
population (Table 6). In accordance with previous stud- 
ies, all individuals were found infected with Wolbachia 
and all of them harbored at least the wCerl strain (Rie- 
gler and Stauffer 2002; Arthofer et al. 2009b, 2011). The 
majority of the Greek individuals were found double- 
infected with wCerl and wCer4 strains while a few 
individuals were found triple-infected with different com- 
binations of wCerl, wCer 2, wCer4, and wCer5 strains or 
infected with four Wolbachia strains (wCerl, wCer2, 
wCer4, and wCer5). It has to be noted, however, that 
tvCerl single-infected individuals were detected in several 
populations while one of them, Lesvos-GR, was found to 
be 100% single-infected with this strain. The majority of 
Krasnodar-RUS individuals were double-infected with 
wCerl and wCer4 strains while triple-infected (wCerl, 
wCer2 and wCer4) flies were also present. The German 
R. cerasi populations were almost fixed for the presence 
of four different Wolbachia strains: wCerl, wCer2, wCer4, 
and wCer5. None of the individuals tested was found 
infected with the wCer3 strain. 

Discussion 

Rhagoletis cerasi present some unique properties when 
compared with the medfiy and the olive fly. Medfly is a 
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Table 6. PCR screening of Rhagoletis cerasi populations with the primer pairs developed for the detection of Wolbachia wCer 1-5 strains (Artho- 
fer et al. 2009b). 
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Wolbachia's presence was screened with the universal Wolbachia 16s rRNA primer pairs (wspec F - 
PCR, polymerase chain reaction; N, number of flies screened. 

Grey shading highlights populations from Germany (light gray) and Russia (dark gray). 



wspecR). 



multivoltine and polyphagous species (Papadopoulos 
et al. 2001), the olive fly is a multivoltine and monopha- 
gous species (Tzanakakis 2003), while the cherry fruit fly 
is a univoltine and stenophagous species (Moraiti et al. 
2012a,b). These three species also differ in their Wolbachi- 
a infection status: natural populations of olive fly and 
medfly are Wolbachia free (Bourtzis et al. 1994; Zabalou 
et al. 2004; Apostolaki et al. 2011; but see Rocha et al. 
2005 for an infected Brazilian population) while Wolba- 
chia has been detected in all natural populations of the 
European cherry fruit fly tested so far. Our results indi- 
cate the presence of structuring in natural populations of 
R. cerasi, which, especially for Greece, is not in accor- 
dance with data for the olive fly (Augustinos et al. 2005; 
Zygouridis et al. 2009), and the medfly (K. Economou 
et al. unpubl. data), which point to rather homogeneous 
groups for both. Is then the genetic structuring observed 
in natural cherry fruit fly populations due to ecological 
factors or rather due to the presence of different/multiple 
Wolbachia strains? 

Variability of microsatellite markers and 
HWE 

We used 13 microsatellite markers; 10 cross-amplified 
from other Rhagoletis species (see Augustinos et al. 
2011a), one from the olive fly (Augustinos et al. 2008), 
and two of the set that has been developed de novo for 



R. cerasi (Arthofer et al. 2009a). Some of the cross-ampli- 
fied markers identified low polymorphism compared with 
species in which they were developed. This is not surpris- 
ing as it is known that cross-amplification of microsatel- 
lite markers includes three major risks: (a) amplification 
of nonhomologous regions, (b) amplification of homolo- 
gous regions with fewer or disrupted repeats, leading 
therefore to lower mutation rates, and (c) poor amplifica- 
tion in some or most of the individuals, due to changes 
in primer binding sites. Sequencing analysis of cross- 
amplified products verified that we amplified homologous 
regions (data not shown), with the exception of RcMic 
(Boms3b) from the olive fly and RcRpl from R. pomonel- 
la. In most of the cases, microsatellites with fewer and/or 
interrupted repeats or with new motifs (data not shown) 
were present, explaining the reduced polymorphism. 
There was also evidence of size variation not attributed 
to motif changes, even within R. cerasi (data not shown), 
implying the participation of other events, such as 
unequal crossing over and/or indels affecting the micro- 
satellite marker size. Although 13 de novo developed 
microsatellite markers were available for R. cerasi (Artho- 
fer et al. 2009a), we managed to use only two of them. 
In preliminary trials, we found that scoring of the 
remaining 11 markers was very difficult and unreliable, 
also due to secondary bands and poor quality of main 
bands (faint bands or not sharp enough to be scored). 
This is probably because these markers harbor multiple 
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motifs, one of which being one base pair repeat motif, 
making their analysis problematic, especially in polyacryl- 
amide gels. 

Deviations from HWE were found in few of the popu- 
lations for some of the cross-amplified loci but not for 
the two species-specific microsatellite markers (Table 2). 
This could partially be attributed to the small sample size 
(in some cases) or the presence of null alleles, as men- 
tioned above. Evidence for a relative extended presence of 
null alleles in some of the samples was found only for 
Rcmic(Rpll). 

Rhagoletis cerasi natural populations' 
polymorphism 

All populations exhibited a comparable degree of poly- 
morphism. The most polymorphic population is the one 
from Krasnodar-RUS. This can be explained considering 
that R. cerasi moved to Europe from western Asia through 
the Caucasus. This population should therefore be closest 
to the original expansion area of the species. Populations 
near the expansion center of a species are expected to be 
more polymorphic, if no drastic bottlenecks have occurred 
and populations have maintained a relative large size 
through time. On the other hand, populations from the 
Eastern Aegean islands (Chios and Lesvos) along with 
Kernitsa from Peloponnesus present the lowest polymor- 
phism. This can be attributed to an island effect, as 
isolated populations are expected to be less polymorphic, 
due to restriction of gene flow and probably small effective 
population size. As far as the sample from Kernitsa (Pelo- 
ponnesus) is concerned, Peloponnesus is not a true island 
and communicates with the rest of mainland Greece; how- 
ever, cherry orchards in the area are few and dispersed. 
Another factor to be taken into account is that the sample 
from Kernitsa consists of individuals that have undergone 
prolonged dormancy, while all other samples consist of 
flies that emerged after annual dormancy. Although the 
prolonged dormancy in R. cerasi pupae seems to be a plas- 
tic response to environmental cues (C. A. Moraiti and N. 
T. Papadopoulos unpubl. data), its possible impact on 
population structuring including bottleneck phenomena 
deserves further investigation. 

Structuring of populations 

Despite the low overall polymorphism observed, the 
markers used here identified well and poorly separated 
clusters (Figs. 2, 3, 4, 5). The well separated groups 
included populations from Aegean islands-GR and 
Krasnodar-RUS followed by Kernitsa-GR and German 
samples. The two loose groups were those from Macedo- 
nia/Epirus-GR and Thessaly-GR. Despite its geographical 
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isolation, the single population from Chania was not well 
differentiated probably due to trade and airport connec- 
tions between Crete and many mainland regions in 
Greece and other European countries including Germany. 

The genetic difference of the samples collected from 
islands or "island"-like sites (Lesvos, Chios, Kernitsa) is 
accompanied by a slightly lower degree of polymorphism 
(Table 1). This is in accordance with the hypothesis of 
partially isolated populations. No evidence of recent bot- 
tlenecks was found (with the exception of Kernitsa popu- 
lation discussed above) suggesting that we are probably 
dealing with long-standing populations that manage to 
keep a relative high population density through time. 

Another finding, although with low support (Figs. 3, 
5), is the topology of the Macedonian-GR group of sam- 
ples (especially the sample from Kastoria): as said before, 
they form a loose group, placed near the Thessaly-GR 
cluster, having, however, a tendency to link with the 
German cluster. If we take into account that Thessaly is 
in Central Greece and Macedonia-GR is mainly in Northern 
Greece and is the region that neighbors with Albania, 
Former Yugoslav Republic of Macedonia and Bulgaria, we 
can speculate in favor of the existence of gene flow from 
Northern Europe in this area. The lack of samples from 
this area limits our analysis. 

Samples from the two Eastern Aegean islands are 
strongly differentiated from all other Greek and German 
samples. They seem to form a "loose" cluster with the 
sample from Russia, pointing toward gene flow from 
Asia Minor, most probably through Turkey, a major 
cherry producer. Analysis of samples from the coastal 
area of Turkey all the way to the coast of Black Sea 
would definitely shed light on gene flow from Asian to 
European R. cerasi populations. It is worth noting that 
these island populations do not form one panmictic pop- 
ulation, despite the fact that the distance between them 
is only 90 km. 

Despite differences in genetic distance within the same 
geographic scale between R. cerasi, C. capitata, and 
B. oleae, the degree of polymorphism was comparable or 
even smaller in R. cerasi (though the same markers 
couldn't be used). Therefore, these differences are not 
marker biased. They are also not population biased, as rel- 
atively large samples were used in all studies (20-30 indi- 
viduals per population). It seems that the fragmented 
landscape of mainland Greece, together with a patchy dis- 
tribution of cherry orchards, the low dispersion ability of 
the cherry fly and local adaptation from close synchroni- 
zation of adult flight with the sweet cherry ripening per- 
iod may restrict gene flow to relatively short distances 
(Phillips and Dirks 1933; Jones and Wallace 1955; Boiler 
and Prokopy 1976; Fletcher 1989; Kneifl et al. 1997). 
Indeed, recent evidence demonstrates high differences in 
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diapause intensity between coastal and highland _R. cerasi 
populations (Papanastasiou et al. 2011) and extended dif- 
ference in adult life history traits between populations 
(Moraiti et al. 2012b). Low reproduction rates of R. cerasi 
because of univoltinism and stenophagy may also partially 
explain differences in genetic structure against both 
C. capitata and B. oleae. On the other hand, infection of 
all natural populations of _R. cerasi with different strains 
of Wolbachia may be yet another factor contributing 
toward structuring and genetic isolation. 

Wolbachia infections in R. cerasi 

Wolbachia is assumed to invade natural populations rap- 
idly and is implicated in speciation phenomena and 
restriction of gene flow between natural populations 
(Shoemaker et al. 1999; Bordenstein et al. 2001; Jaenike 
et al. 2006; Koukou et al. 2006; Miller et al. 2010; Branca 
et al. 2011). This can be either through CI and/or rapid 
fixation of specific genotypes during invasion. Natural 
European _R. cerasi populations are infected with multiple 
Wolbachia strains and show different infection patterns. 
For example, in Poland, Italy, and Austria, there are pop- 
ulations infected with all five wCer strains, in the Czech 
Republic and Portugal there are populations infected with 
four of the five strains (except wCer2), while in Switzer- 
land there was an infection with all wCer strains except 
wCer3 (Arthofer et al. 2009b, 2011). Moreover, it has 
been shown that the presence of different Wolbachia 
strains in R. cerasi could be the causal factor (Riegler and 
Stauffer 2002) of the incompatibility observed among 
R. cerasi natural populations in previous studies (Boiler 
et al. 1976; Matolin 1976). 

All individuals in this study (10-20 per population) 
were infected with Wolbachia. The analysis of six individ- 
uals from different samples (genetically or geographically 
distant) through cloning and sequencing analysis, as well 
as the PCR screening of more individuals with the diag- 
nostic PCR primers, pointed toward (a) the presence of 
complex patterns of infections with the four of the five 
known wCer strains (1, 2, 4, and 5) and (b) the possible 
existence of new, still uncharacterized strains (see Chania 
population). However, not all genes resulted into the 
same infection pattern (Tables 5, 6). This can be attrib- 
uted to the preferential recovery of specific sequences 
after PCR and/or cloning. Overall, these two approaches 
(PCR screening and cloning plus sequencing) produced 
similar results and supported the presence of multiple 
infections. However, both PCR and cloning can lead in 
bias in favor of specific sequences and to the misinterpre- 
tation of the "actual" infection status. On the other hand, 
PCR screening allows the analysis at population level but 
currently is limited by the development of strain-specific 
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primer pairs only for the wsp gene and by the inability to 
detect "hidden" strains. 

Arthofer et al. (2011) highlighted the difficulties in 
Wolbachia strain characterization through MLST in multi- 
ple-infected hosts. They also provided a tool to analyze 
complex infection types (Allelic Intersection Analysis) and 
evaluated this tool for wCer Wolbachia strains of R. cerasi 
by assigning specific MLST alleles to four of the five wCer 
strains (Table 5). It is evident that the MLST- and wsp- 
based typing approaches do not have the power to resolve 
complex infection patterns such as those observed in 
R. cerasi. There is probably the need for new tools and 
approaches to address such questions. Current use and 
evaluation of different types of molecular markers, includ- 
ing Variable Number Tandem Repeat markers (VNTRs) 
(Riegler et al. 2005, 2012; Schneider et al. 2013) can 
probably contribute toward this direction. 

Multiple and/or low titer strains have been shown to 
contribute to genetic differentiation of their hosts (Artho- 
fer et al. 2009b, 2011; Schneider et al. 2013). Our data 
confirm the presence of multiple infections, different 
types of infections in different samples and polymorphic 
infection status within the samples. They also raise the 
possibility of the presence of still unknown strains. Given 
also the ability of at least two strains (wCer2 and wCer4) 
to induce CI (Boiler et al. 1976; Riegler et al. 2004; Zaba- 
lou et al. 2004, 2009; Apostolaki et al. 2011), Wolbachia is 
not a factor that can be overlooked when addressing the 
question of population structuring in _R. cerasi. 

Recombination has perhaps been the most significant 
contributing factor to the evolution of Wolbachia and may, 
at least partly, explain the diversity of the wCer strains pres- 
ent in European populations of R. cerasi (Jiggins 2002; 
Baldo and Werren 2007; Arthofer et al. 2009b; Klasson 
et al. 2009). Low titer infections are an additional, hidden 
source of Wolbachia diversity, as documented in different 
hosts (Arthofer et al. 2009b; Miller et al. 2010; Augustinos 
et al. 2011b). There is also a number of studies showing that 
the transfer of Wolbachia genes (in some cases even the 
entire genome) into host chromosomes might be a rather 
common phenomenon (Kondo et al. 2002; Hotopp et al. 
2007; Nikoh et al. 2008; Aikawa et al. 2009; Klasson et al. 
2009; Nikoh and Nakabachi 2009; Woolfit et al. 2009; Dou- 
doumis et al. 2012). If this is the case, then nucleus- 
encoded Wolbachia genes may also contribute to the 
observed genetic diversity. However, our data do not pro- 
vide evidence for either (intragene) recombination or hori- 
zontal gene transfer events. 

Do all wCer Wolbachia strains induce CI? The fact that 
R. cerasi is a univoltine species and an efficient and 
robust lab rearing system is lacking makes any experimen- 
tation to address questions regarding the CI properties of 
wCer strains very difficult. Interpretation of historical 
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data in European populations of R. cerasi as well as trans- 
infection experiments with medfly and olivefly suggest 
that wCer2 and wCer4 induce strong CI into their hosts 
and are bidirectionally incompatible (Boiler et al. 1976; 
Riegler and Stauffer 2002; Riegler et al. 2004; Zabalou 
et al. 2004, 2009; Apostolaki et al. 2011). On the other 
hand, Drosophila simulans was not shown to be a suitable 
host for the establishment and/or the expression of CI for 
either wCerl or wCer2 strains (Riegler et al. 2004). Thus, 
transinfection experiments in suitable hosts like medfly 
can probably unravel the CI properties of wCer3, wCer5, 
and the two potentially new wCer strains. 

Concluding Remarks 

Our results reveal some structuring of natural European 
cherry fruit fly populations, which can probably be gener- 
alized to its whole geographic range. A number of factors 
may have contributed to this structuring, such as habitat 
fragmentation as well as Wolbachia infection. Microsatel- 
lite analysis of more populations obtained from geograph- 
ically distant areas, along with further genetic 
characterization of Wolbachia strains and their CI proper- 
ties are required to understand the population dynamics 
of this economically important species and its associated 
Wolbachia strains more completely. This knowledge will 
be needed before any novel population control method is 
applied on a large scale for this important pest species. 

Author Contribution Box 

Conceived and designed the experiments: K. B., N. T. P., 
A. A. A. Performed the experiments: A. A. A., A. A. K., 
C. A. M. Analyzed the data: K. B., A. A. A., N. T. P., P. 
M.-T. Contributed reagents/materials/analysis tools: K. B., 
N. T. P., P. M. T. Wrote the manuscript: K. B., A. A. A, 
N. T. P., P. M. T. 

Acknowledgments 

We thank Antigone Zacharopoulou (Department of Biol- 
ogy, University of Patras, Greece) for her generous sup- 
port and encouragement throughout this study. Thanks 
are also extended to Ludmilla Vassilieva, Kirsten Koep- 
pler, Byron Katsoyannos, Ilias Kounatidis, Alex Diaman- 
tidis, Stella Papanastasiou, and Charalampos Ioannou for 
either submitting R. cerasi samples or participating in 
collecting the different populations. We also thank Wolf- 
gang Miller and Daniela Schneider for providing control 
DNA samples used in the wCer strain-specific PCR 
screen presented in Table 6. This research was supported 
in part by the European Community's Seventh Frame- 
work Programme CSA-SA_REGPROT-2007-1 (under 



grant agreement no 203590 (K. B.). This study was par- 
tially supported by IKYDA and the EPEAEK II - Pytha- 
goras (cofounded by the European Social and National 
Resources and the Greek ministry of Education) awarded 
to N. T. P. P. M., N. T. P., and K. B. benefitted from 
travel grants in the frame of the EU COST Action 
FA0701 "Arthropod Symbiosis: from fundamental studies 
to pest and disease management". The authors are grate- 
ful to the FAO/IAEA Coordinated Research Program 
"Use of symbiotic bacteria to reduce mass-rearing costs 
and increase mating success in selected fruit pests in sup- 
port of SIT application" for the overall support of this 
study. The authors also thank Stefan Oehler for his com- 
ments on an earlier version of the manuscript. The fund- 
ing agencies had no role in study design, data collection 
or analysis, decision to publish, or preparation of the 
manuscript. 

Conflict of Interest 

None declared. 
References 

Aikawa, T., H. Anbutsu, N. Nikoh, T. Kikuchi, F. Shibata, and 
T. Fukatsu. 2009. Longicorn beetle that vectors pinewood 
nematode carries many Wolbachia genes on an autosome. 
Proc. Biol. Sci. 276:3791-3798. 

Aketarawong, N., M. Bonizzoni, S. Thanaphum, L. M. 
Gomulski, G. Gasperi, A. R. Malacrida, et al. 2007. 
Inferences on the population structure and colonization 
process of the invasive oriental fruit fly, Bactrocera dorsalis 
(Hendel). Mol. Ecol. 16:3522-3532. 

Apostolaki, A., I. Livadaras, A. Saridaki, A. Chrysargyris, C. 
Savakis, and K. Bourtzis. 2011. Transinfection of the olive 
fruit fly Bactrocera oleae with Wolbachia: towards a 
symbiont-based population control strategy. J. Appl. 
Entomol. 13:546-553. 

Arthofer, W., S. Krumboeck, H. Schuler, B. Rasool, M. Riegler, 
K. Koeppler, et al. 2009a. Thirteen new microsatellite loci in 
Rhagoletis cerasi (Diptera: Tephritidae), a model host species 
for Wolbachia symbiosis in field populations. Mol. Ecol. 
Resour. 9:1460-1466. 

Arthofer, W., M. Riegler, D. Schneider, M. Krammer, 
W. J. Miller, and C. Stauffer. 2009b. Hidden Wolbachia 
diversity in field populations of the European cherry fruit fly, 
Rhagoletis cerasi (Diptera, Tephritidae). Mol. Ecol. 18: 
3816-3830. 

Arthofer, W., M. Riegler, H. Schuler, D. Schneider, K. Moder, 
W. J. Miller, et al. 2011. Allele intersection analysis: a novel 
tool for multi locus sequence assignment in multiply 
infected hosts. PLoS One 6:e22198. 

Atyame, C. M., N. Pasteur, E. Dumas, P. Tortosa, 
M. L. Tantely, N. Pocquet, et al. 2011. Cytoplasmic 



1958 



© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



A. A. Augustinos ef al. 



R. cerasi Microsatellite and Wolbachia Analysis 



incompatibility as a means of controlling Culex pipiens 
quinquefasciatus mosquito in the islands of the south- 
western Indian Ocean. PLoS Negl. Trop. Dis. 5:el440. 

Augustinos, A. A., Z. Mamuris, E. Stratikopoulos, S. D'Amelio, 
A. Zacharopoulou, and K. D. Mathiopoulos. 2005. 
Microsatellite analysis of olive fly populations in the 
Mediterranean indicates a westward expansion of the 
species. Genetica 125:231-241. 

Augustinos, A. A., E. E. Stratikopoulos, E. Drosopoulou, 
E. G. Kakani, P. Mavragani-Tsipidou, A. Zacharopoulou, 
et al. 2008. Isolation and characterization of microsatellite 
markers from the olive fly, Bactrocera oleae, and their cross- 
species amplification in the Tephritidae family. BMC 
Genomics 9:e618. 

Augustinos, A. A., A. K. Asimakopoulou, N. T. Papadopoulos, 
and K. Bourtzis. 201 la. Cross-amplified microsatellites in the 
European cherry fly, Rhagoletis cerasi: medium polymorphic- 
highly informative markers. Bull. Entomol. Res. 101:45-52. 

Augustinos, A. A., D. Santos-Garcia, E. Dionyssopoulou, and 
K. Bourtzis. 2011b. Detection and characterization of 
Wolbachia infections in natural populations of aphids: is the 
hidden diversity fully unravelled? PLoS One 6:e28695. 

Baldo, L., and J. H. Werren. 2007. Revisiting Wolbachia 
supergroup typing based on wsp: spurious lineages and 
discordance with MLST. Curr. Microbiol. 55:81-87. 

Baldo, L., J. C. Dunning Hotopp, K. A. Jolley, S. R. 
Bordenstein, S. A. Biber, R. R. Choudhury, et al. 2006. 
Multilocus sequence typing system for the endosymbiont 
Wolbachia pipientis. Appl. Environ. Microbiol. 72:7098-7110. 

Boiler, E. F., and B. I. Katsoyannos. 1977. Measuring mating 
frequency in Rhagoletis cerasi. Bull. SROP Int. Organ Biol. 
Control Noxious Anim. Plants West Palearctic Reg. Sec. 
5:102-103. 

Boiler, E. F., and R. J. Prokopy. 1976. Bionomics and 

management of Rhagoletis. Annu. Rev. Entomol. 21:223-246. 

Boiler, E. F., A. Haisch, and R. J. Prokopy. 1971. Sterile insect 
release method against Rhagoletis cerasi L.: preparatory 
ecological and behavioural studies. P. 77 Proceedings of 
Symposium Athens, 1970. Sterile principles for insect control 
or eradication. IAEA, Vienna. 

Boiler, E. F., K. Russ, V. Vallo, and G. L. Bush. 1976. 

Incompatible races of European cherry fruit fly, Rhagoletis 
cerasi (Diptera: Tephritidae), their origin and potential use 
in biological control. Entomol. Exp. Appl. 20:237-247. 

Boiler, E. F., U. Remund, and B. I. Katsoyannos. 1977. 
Measuring mating propensity (speed) in Rhagoletis cerasi. 
Pp. 96-97 in E. F.Boller and D. L. Chambers, eds. Quality 
control. An ideal book for fruit fly workers. International 
Organization for Biological Control of Noxious Animals and 
Plants, WPRS Bulletin 1977/5. 

Boiler, E. F., B. I. Katsoyannos, and C. Hippe. 1998. Host 
races of Rhagoletis cerasi L. (Diptera: Tephritidae): effect of 
prior adult experience on oviposition site preference. 
J. Appl. Entomol. 122:231-237. 



2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



Bonizzoni, M., B. I. Katsoyannos, R. Marguerie, 

C. R. Guglielmino, G. Gasperi, A. Malacrida, et al. 2002. 
Microsatellite analysis reveals remating by wild 
Mediterranean fruit fly females, Ceratitis capitata. Mol. Ecol. 
11:1915-1921. 

Bonizzoni, M., C. R. Guglielmino, C. J. Smallridge, 

M. Gomulski, A. R. Malacrida, and G. Gasperi. 2004. On 
the origins of medfly invasion and expansion in Australia. 
Mol. Ecol. 13:3845-3855. 

Bordenstein, S. R., F. P. O'Hara, and J. H. Werren. 2001. 
Wolbachia-induced incompatibility precedes other hybrid 
incompatibilities in Nasonia. Nature 409:707-710. 

Bourtzis, K., and A. S. Robinson. 2006. Insect pest control 
using Wolbachia and/or radiation. Pp. 225-246 in K. 
Bourtzis, T. Miller, eds. Insect symbiosis 2. CRC Press, 
Talylor and Francis Group, LLC, FL. 

Bourtzis, K., A. Nirgianaki, P. Onyango, and C. Savakis. 1994. 
A prokaryotic dnaA sequence in Drosophila melanogaster: 
Wolbachia infection and cytoplasmic incompatibility among 
laboratory strains. Insect Mol. Biol. 3:131-142. 

Branca, A., B. P. Le Ru, F. Vavre, J.-F. Silvain, and S. Dupas. 
2011. Intraspecific specialization of the generalist parasitoid 
Cotesia sesamiae revealed by polyDNAvirus polymorphism 
and associated with different Wolbachia infection. Mol. Ecol. 
20:959-971. 

Brelsfoard, C. L., and S. L. Dobson. 2009. Wolbachia-based 
strategies to control insect pests and disease vectors. Asia 
Pac. J. Mol. Biol. Biotechnol. 17:55-63. 

Casiraghi, M., S. R. Bordenstein, L. Baldo, N. Lo, T. Beninati, 
J. J. Wernegreen, et al. 2005. Phylogeny of Wolbachia 
pipientis based on gltA, groEL and ftsZ gene sequences: 
clustering of arthropod and nematode symbionts in the F 
supergroup, and evidence for further diversity in the 
Wolbachia tree. Microbiology 151:4015^022. 

Chen, Y. H., S. H. Berlocher, S. B. Opp, and G. K. Roderick. 
2010. Post-colonization temporal genetic variation of an 
introduced fly, Rhagoletis completa. Genetica 138:1059-1075. 

Cornuet, J. M., and G. Luikart. 1996. Description and power 
analysis of two tests for detecting recent population 
bottlenecks from allele frequency data. Genetics 
144:2001-2014. 

Daniel, C, and J. Grunder. 2012. Inergrated management of 
European cherru fruit fly Rhagoletis cerasi (L.): situation in 
Switzerland and Europe. Insects 3:956-988. 

Doudoumis, V., G. Tsiamis, F. Wamwiri, C. Brelsfoard, 
U. Alam, and E. Aksoy. 2012. Detection and 
characterization of Wolbachia infections in laboratory and 
natural populations of different species of tsetse flies (genus 
Glossina). BMC Microbiol. 12(Suppl. l):eS3. 

Drew, R. A. I., J. Ma, S. Smith, and J. M. Hughes. 2011. The 
taxonomy and phylogenetic relationships of species in the 
Bactrocera musae complex of fruit flies (Diptera: 
Tephritidae: Dacinae) in Papua New Guinea. Raffles Bull. 
Zool. 59:145-162. 



1959 



R. cerasi Microsatellite and Wolbachia Analysis 

Drosopoulou, E., K. Koeppler, I. Kounatidis, I. Nakou, 
N. T. Papadopoulos, K. Bourtzis, et al. 2010. Genetic and 
cytogenetic analysis of the Walnut-Husk fly (Diptera: 
Tephritidae). Ann. Entomol. Soc. Am. 103:1003-1011. 

Drosopoulou, E., A. A. Augustinos, I. Nakou, K. Koeppler, 
I. Kounatidis, H. Vogt, et al. 2011. Genetic and cytogenetic 
analysis of the American cherry fruit fly, Rhagoletis cingulata 
(Diptera: Tephritidae). Genetica 139:1449-1464. 

Evanno, G., S. Regnaut, and J. Goudet. 2005. Detecting the 
number of clusters of individuals using the software 
STRUCTURE: a simulation study. Mol. Ecol. 14: 
2611-2620. 

Felsenstein, J. 1994. PHYLIP (Phylogeny Inference Package) 
version 3.6. Distributed by the Author. Department of 
Genome Sciences, University of Washington, Seattle. 

Fimiani, P. 1989. Mediterranean region. Pp. 39-50 in A. S. 
Robinson and G. Hopper, eds. Fruit flies: their biology, 
natural enemies and control. Elsevier, Amsterdam, 
Netherlands. 

Fletcher, B. S. 1989. Life history strategies of tephritid flies. Pp. 

209-219 in A. S. Robinson and G. Hooper, eds. Fruit flies: 

their biology, natural enemies and control World crop pests 

3(B). Elsevier, Amsterdam, Netherlands. 
Forbes, A. A., T. H. Q. Powell, L. L. Stelinski, J. J. Smith, and 

J. L. Feder. 2009. Sequential sympatric speciation across 

trophic levels. Science 323:776-779. 
Goudet, J. 2001. FSTAT, a program to estimate and test gene 

diversities and fixation indices (version 2.9.3). 
Guo, S. W., and E. A. Thomson. 1992. Performing the exact 

test of Hardy-Weinberg proportions for multiple alleles. 

Biometrics 48:371-372. 
Hotopp, J. C., M. E. Clark, D. C. Oliveira, J. M. Foster, 

P. Fischer, and M. C. Muhoz Torres. 2007. Widespread 

lateral gene transfer from intracellular bacteria to 

multicellular eukaryotes. Science 317:1753-1756. 
Jaenike, J., K. A. Dyer, C. Cornish, and M. S. Minhas. 2006. 

Asymmetrical reinforcement and Wolbachia infection in 

Dwsophila. PLoS Biol. 4:1852-1862. 
Jiggins, F. M. 2002. The rate of recombination in Wolbachia 

bacteria. Mol. Biol. Evol. 19:1640-1643. 
Jones, S. C, and L. Wallace. 1955. Cherry fruit fly dispersion 

studies. J. Econ. Entomol. 48:616-617. 
Katsoyannos, B. I., N. T. Papadopoulos, and D. Stavridis. 

2000. Evaluation of trap types and food attractants for 

Rhagoletis cerasi (Diptera: Tephritidae). J. Econ. Entomol. 

93:1005-1010. 

Khamis, F. M., N. Karam, S. Ekesi, M. DE Meyer, A. Bonomi, 
and L. M. Gomulski. 2009. Uncovering the tracks of a 
recent and rapid invasion: the case of the fruit fly pest 
Bactrocera invadens (Diptera: Tephritidae) in Africa. Mol. 
Ecol. 18:4798-4810. 

Klasson, L., Z. Kambris, E. P. Cook, T. Walker, and 

S. P. Sinkins. 2009. Horizontal gene transfer between Wolbachia 
and the mosquito Aedes aegypti. BMC Genomics 10:33. 



A. A. Augustinos ef al. 

Kneifl, V., F. Paprstein, and J. Khourkova. 1997. Dispersion of 
cherry fruit fly (Rhagoletis cerasi L.). Vedecke Prace 
Ovocnarsk 15:89-92. 

Kondo, N, N. Nikoh, N. Ijichi, M. Shimada, and T. Fukatsu. 
2002. Genome fragment of Wolbachia endosymbiont 
transferred to X chromosome of host insect. PNAS 
99:14280-14285. 

Koukou, K., H. Pavlikaki, G. Kilias, J. H. Werren, 
K. Bourtzis, and S. N. Alahiotis. 2006. Influence of 
antibiotic treatment and Wolbachia curing on sexual 
isolation among Drosophila melanogaster cage populations. 
Evolution 60:87-96. 

Kounatidis, I., N. Papadopoulos, K. Bourtzis, and 

P. Mavragani-Tsipidou. 2008. Genetic and cytogenetic 
analysis of the fruit fly Rhagoletis cerasi (Diptera: 
Tephritidae). Genome 51:479-491. 

Kovanci, O. B., and B. Kovanci. 2006. Effect of altitude on 
seasonal flight activity of Rhagoletis cerasi flies (Diptera: 
Tephritidae). Bull. Entomol. Res. 96:345-351. 

Martin, D. P., P. Lemey, M. Lott, V. Moulton, D. Posada, and 
P. Lefeuvre. 2010. RDP3: a flexible and fast computer 
program for analyzing recombination. Bioinformatics 
26:2462-2463. 

Matolin, S. 1976. Mechanism causing incompatibility between 
different strains of Rhagoletis cerasi (Diptera, Tephritidae). 
Acta Entomol. Bohemoslov. 73:289-292. 

Maxwell, S. A., G. Rasic, and N. Keyghobadi. 2009. 

Characterization of microsatellite loci for the western cherry 
fruit fly, Rhagoletis indifferens (Diptera: Tephritidae). Mol. 
Ecol. Resour. 9:1025-1028. 

Meixner, M. D., B. A. Mcpheron, J. G. Silva, G. E. Gasparich, 
and W. S. Sheppard. 2002. The Mediterranean fruit fly in 
California: evidence for multiple introductions and 
persistent populations based on microsatellite and 
mitochondrial DNA variability. Mol. Ecol. 11:891-899. 

Michel, A. P., J. Rull, M. Aluja, and J. L. Feder. 2007. The 
genetic structure of hawthorn-infesting Rhagoletis pomonella 
populations in Mexico: implications for sympatric host race 
formation. Mol. Ecol. 16:2867-2878. 

Michel, A. P., S. Sim, T. H. Q. Powell, M. S. Taylor, P. Nosil, 
and J. L. Feder. 2010. Widespread genomic divergence 
during sympatric speciation. PNAS 107:9724-9729. 

Miller, W. J., L. Ehrman, and D. Schneider. 2010. Infectious 
speciation revisited: impact of symbiont-depletion on female 
fitness and mating behavior of Drosophila paulistoru. PLoS 
Pathog. 6:el001214. 

Moraiti, C. A., C. T. Nakas, and N. T. Papadopoulos. 2012a. 
Prolonged pupal dormancy is associated with significant 
fitness cost for adults of Rhagoletis cerasi (Diptera: 
Tephritidae). J. Insect Physiol. 58:1128-1135. 

Moraiti, C. A., C. T. Nakas, K. Koppler, and 

N. T. Papadopoulos. 2012b. Geographical variation in adult 
life history traits of the European cherry fruit fly, Rhagoletis 
cerasi (Diptera: Tephritidae). Biol. J. Linn. Soc. 107:137-152. 



1960 



© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



A. A. Augustinos ef al. 



R. cerasi Microsatellite and Wolbachia Analysis 



Nardi, F., A. Carapelli, R. Dallai, G. K. Roderick, and F. Frati. 
2005. Population structure and colonization history of the 
olive fly, Bactrocera oleae (Diptera, Tephritidae) . Mol. Ecol. 
14:2729-2738. 

Nei, M. 1972. Genetic distance between populations. Am. Nat. 
106:283-292. 

Nikoh, N., and A. Nakabachi. 2009. Aphids acquired symbiotic 
genes via lateral gene transfer. BMC Biol. 7:el2. 

Nikoh, N., K. Tanaka, F. Shibata, N. Kondo, M. Hizume, M. 
Shimada, et al. 2008. Wolbachia genome integrated in an 
insect chromosome: evolution and fate of laterally 
transferred endosymbiont genes. Genome Res. 18:272-280. 

Papadopoulos, N. T., B. I. Katsoyannos, J. R. Carey, and 
N. A. Kouloussis. 2001. Seasonal and annual occurrence of 
the Mediterranean fruit fly (Diptera: Tephritidae) in 
Northern Greece. Ann. Entomol. Soc. Am. 94:41-50. 

Papanastasiou, S. A., D. Nestel, A. D. Diamantidis, C. T. 
Nakas, and N. T. Papadopoulos. 2011. Physiological and 
biological patterns of a highland and a coastal population of 
the European cherry fruit fly during diapause. J. Insect 
Physiol. 57:83-93. 

Paraskevopoulos, C, S. R. Bordenstein, J. J. Wernegreen, 
J. H. Werren, and K. Bourtzis. 2006. Toward a Wolbachia 
multilocus sequence typing system: discrimination of 
Wolbachia strains present in Drosophila species. Curr. 
Microbiol. 53:388-395. 

Peakall, R., and P. E. Smouse. 2006. GENALEX 6: genetic 
analysis in Excel. Population genetic software for teaching 
and research. Mol. Ecol. Notes 6:288-295. 

Phillips, C. R., and C. O. Dirks. 1933. Dispersal of the apple 
maggot. J. Econ. Entomol. 26:344-349. 

Piry, S., A. Alapetite, J. Cornuet, D. Paetkau, L. Baudouin, and 
A. Estoup. 2004. GENECLASS2: a software for genetic 
assignment and first-generation migrant detection. J. Hered. 
95:536-539. 

Pritchard, J. K., M. Stephens, and P. Donnelly. 2000. Inference 

of population structure using multilocus genotype data. 

Genetics 155:945-959. 
Rannala, B., and J. L. Mountain. 1997. Detecting immigration 

by using multilocus genotypes. PNAS 94:9197-9201. 
Raymond, M., and F. Rousset. 1995. GENEPOP v. 1.2: 

population genetics software for exact tests and 

ecumenicism. J. Hered. 86:248-249. 
Riegler, M., and C. Stauffer. 2002. Wolbachia infections and 

superinfections in cytoplasmically incompatible populations 

of the European cherry fruit fly Rhagoletis cerasi (Diptera, 

Tephritidae). Mol. Ecol. 11:2425-2434. 
Riegler, M., S. Charlat, C. Stauffer, and H. Mercot. 2004. 

Wolbachia transfer from Rhagoletis cerasi to Drosophila 

simulans: investigating the outcomes of host-symbiont 

coevolution. Appl. Environ. Microbiol. 70:273-279. 
Riegler, M., M. Sidhu, W. J. Miller, and S. L. O'Neil. 2005. 

Evidence for a global Wolbachia replacement in Drosophila 

melanogaster. Curr. Biol. 15:1428-1433. 



Riegler, M., I. Iturbe-Ormaetxe, M. Woolfit, W. J. Miller, and 
S. L. O'Neil. 2012. Tandem repeat markers as novel 
diagnostic tools for high resolution fingerprinting of 
Wolbachia. BMC Microbiol. 12:S12. 

Rocha, L. S., R. O. Mascarenhas, A. L. P. Perondini, and D. 
Selivon. 2005. Occurrence of Wolbachia in Brazilian samples 
of Ceratitis capitata (Wiedemann) (Diptera: Tephritidae). 
Neotrop. Entomol. 34:1013-1015. 

Ros, V. I., V. M. Fleming, E. J. Feil, and J. A. Breeuwer. 2009. 
How diverse is the genus Wolbachia? Multiple-gene 
sequencing reveals a putatively new Wolbachia supergroup 
recovered from spider mites (Acari: Tetranychidae). Appl. 
Environ. Microbiol. 75:1036-1043. 

Rousset, F., and M. Raymond. 1995. Testing heterozygote 
excess and deficiency. Genetics 140:1413-1419. 

Saridaki, A., and K. Bourtzis. 2010. Wolbachia: more than just 
a bug in insects genitals. Curr. Opin. Microbiol. 13:67-72. 

Schneider, D. I., K. I. Garschall, A. G. Parker, A. M. M. Abd- 
AUa, and W. J. Miller. 2013. Global Wolbachia prevalence, 
titer fluctuations and their potential of causing cytoplasmic 
incompatibilities in tsetse flies and hybrids of Glossina 
morsitans subgroup species. J. Invertebr. Pathol. 112:104- 
115. 

Schuler, H, W. Arthofer, S. Krumbock, K. Koppler, H. Vogt, 
L. A. F. Texeira, et al. 2009. The bacterial endosymbiont 
Wolbachia in the invasive cherry fruit fly Rhagoletis cingulata 
(Diptera, Tephritidae). Mitteilungen der Deutschen Gesellschaft 
fu"r Allgemeine und Angewandte Entomologie 17:99-101. 

Schuler, H, W. Arthofer, M. Riegler, C. Bertheau, S. 
Krumborck, K. Koppler, et al. 2011. Multiple Wolbachia 
infections in Rhagoletis pomonella. Entomol. Exp. Appl. 
139:138-144. 

Schwarz, D., B. A. McPheron, G. B. Hartl, E. F. Boiler, and 
T. S. Hoffmeister. 2003. A second case of genetic host races 
in Rhagoletis? A population genetic comparison of 
sympatric host populations in the European cherry fruit fly, 
Rhagoletis cerasi. Entomol. Exp. Appl. 108:11-17. 

Serbus, L. R., C. Casper-Lindley, F. Landmann, and 
W. Sullivan. 2008. The genetics and cell biology of 
Wolbachia-host interactions. Annu. Rev. Genet. 42:683-707. 

Shoemaker, D. D., V. Katju, and J. Jaenike. 1999. Wolbachia 
and the evolution of reproductive isolation between 
Drosophila recens and Drosophila subquinaria. Evolution 
53:1157-1164. 

Song, S. D., R. A. I. Drew, and J. M. Hughes. 2007. Multiple 
paternity in a natural population of a wild tobacco fly, 
Bactrocera cacuminata (Diptera: Tephritidae), assessed by 
microsatellite DNA markers. Mol. Ecol. 16:2353-2361. 

Tamura, K., D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. 
Kumar. 2011. MEGA5: molecular evolutionary genetics 
analysis using maximum likelihood, evolutionary distance, and 
maximum parsimony methods. Mol. Biol. Evol. 28:2731-2739. 

Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. 
CLUSTAL W: improving the sensitivity of progressive 



2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



1961 



R. cerasi Microsatellite and Wolbachia Analysis 



A. A. Augustinos ef al. 



multiple sequence alignment through sequence weighting, 
position-specific gap penalties and weight matrix choice. 
Nucleic Acids Res. 22:4673-4680. 

Tzanakakis, M. E. 2003. Seasonal development and dormancy 
of insects and mites feeding on olive: a review. Neth. J. 
Zool. 52:87-224. 

Velez, S., M. S. Taylor, M. A. F. Noor, N. F. Lobo, and 
J. L. Feder. 2006. Isolation and characterization of 
microsatellite loci from the apple maggot fly 
Rhagoletis pomonella (Diptera: Tephritidae). Mol. Ecol. 
Notes 6:90-92. 

Virgilio, M., H. Delatte, T. Backeljau, and M. De Meyer. 2010. 

Macrogeographic population structuring in the 

cosmopolitan agricultural pest Bactrocera cucurbitae 

(Diptera: Tephritidae). Mol. Ecol. 19:2713-2724. 
Walker, T., P. H. Johnson, L. A. Moreira, I. Iturbe-Ormaetxe, 

F. D. Frentiu, and C. J. McMeniman. 2011. The wMel 

Wolbachia strain blocks dengue and invades caged Aedes 

aegypti populations. Nature 476:450-455. 
Werren, J. H., and D. M. Windsor. 2000. Wolbachia infection 

frequencies in insects: evidence of a global equilibrium? 

Proc. Biol. Sci. 267:1277-1285. 
Werren, J. H., L. Baldo, and M. E. Clark. 2008. Wolbachia: 

master manipulators of invertebrate biology. Nat. Rev. 

Microbiol. 6:741-751. 
White, I. M., and M. M. Elson-Harris. 1992. Fruit flies of 

economic significance: their identification and bionomics. 

CAB International, Wallinford, UK. 
Woolfit, M., I. Iturbe-Ormaetxe, E. A. McGraw, and 

S. L. O'Neill. 2009. An ancient horizontal gene transfer 

between mosquito and the endosymbiotic bacterium 

Wolbachia pipientis. Mol. Biol. Evol. 26:367-374. 
Yeh, F. C, T. Boyle, Y. Rongcai, Z. Ye, and J. M. Xiyan. 1999. 

POPGENE. 

Zabalou, S., M. Riegler, M. Theodorakopoulou, C. Stauffer, 
C. Savakis, and K. Bourtzis. 2004. Wb/bac/na-induced 
cytoplasmic incompatibility as a means for insect pest 
population control. PNAS 101:15042-15045. 

Zabalou, S., A. Apostolaki, I. Livadaras, G. Franz, 

A. S. Robinson, C. Savakis, et al. 2009. Incompatible insect 
technique: incompatible males from a Ceratitis capitata 
(Diptera: Tephritidae) genetic sexing strain. Entomol. Exp. 
Appl. 132:232-240. 

Zygouridis, N. E., A. A. Augustinos, F. G. Zalom, and 

K. D. Mathiopoulos. 2009. Analysis of olive fly invasion in 
California based on microsatellite markers. Heredity 102: 
402^12. 



Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Figure SI. Posterior probability of possible number of pop- 
ulations in STRUCTURE. Note the potential presence of six 
genetic clusters indicated by the arrow in diagram (D). 
Figure S2. Maximum Likelihood dendrogram based on the 
16S rRNA gene sequences. Gray shade: sequences generated 
in this study. Wolbachia 16S rRNA gene sequences repre- 
senting current known Supergroups were retrieved from 
Genebank. For these sequences, the host scientific name, 
the Genebank Accession Number and the respective super- 
group are shown. Analyses were conducted in Mega 5 
(Tamura et al. 2011). The evolutionary history was inferred 
using the Maximum Likelihood method based on the Tam- 
ura-Nei model. The tree with the highest log likelihood is 
shown and percentages of trees in which the associated taxa 
clustered together are indicated next to the branches. 
Figure S3. Maximum Likelihood dendrogram based on the 
coxA gene sequences. Gray shade: sequences generated in 
this study. Darker gray shade and bold: sequences of known 
wCer strains. Different Wolbachia coxA alleles were 
retrieved from the Wolbachia MLST database. Alleles 
representing all available supergroups were chosen. New 
alleles coxA- 188 and coxA- 189 are marked with asterisks 
(*). Analyses were conducted in Mega 5 (Tamura et al. 
2011). The evolutionary history was inferred using the 
Maximum Likelihood method based on the Tamura-Nei 
model. The tree with the highest log likelihood is shown 
and percentages of trees in which the associated taxa clus- 
tered together are indicated next to the branches. 
Figure S4. Maximum Likelihood dendrogram based on 
hcpA gene sequences. Gray shade: sequences generated in 
this study. Darker gray shade and bold: sequences of known 
wCer strains. Different Wolbachia hcpA alleles were 
retrieved from the Wolbachia MLST database. Alleles repre- 
senting all available supergroups were chosen. New alleles 
hpA-212 and hcpA-213 are marked with asterisks (*). Anal- 
yses were conducted in Mega 5 (Tamura et al. 2011). The 
evolutionary history was inferred using the Maximum Like- 
lihood method based on the Tamura-Nei model. The tree 
with the highest log likelihood is shown and percentages of 
trees in which the associated taxa clustered together are 
shown next to the branches. 
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